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A turbulent flow is maintained by an external supply of kinetic energy, which is eventually dissipated into 
heat at steep velocity gradients. The scale at which energy is supplied greatly differs from the scale at which 
energy is dissipated, the more so as the turbulent intensity (the Reynolds number) is larger. The resulting 
energy flux over the range of scales, intermediate between energy injection and dissipation, acts as a source of 
time irreversibility. As it is now possible to follow accurately fluid particles in a turbulent flow field, both 
from laboratory experiments and from numerical simulations, a natural question arises: how do we detect 
time irreversibility from these Lagrangian data? Here we discuss recent results concerning this problem. For 
Lagrangian statistics involving more than one fluid particle, the distance between fluid particles introduces an 
intrinsic length scale into the problem. The evolution of quantities dependent on the relative motion between 
these fluid particles, including the kinetic energy in the relative motion, or the configuration of an initially 
isotropic structure can be related to the equal-time correlation functions of the velocity field, and is therefore 
sensitive to the energy flux through scales, hence to the irreversibility of the flow. In contrast, for single-particle 
Lagrangian statistics, the most often studied velocity structure functions cannot distinguish the “arrow of time”. 
Recent observations from experimental and numerical simulation data, however, show that the change of kinetic 
energy following the particle motion, is sensitive to time-reversal. We end the survey with a brief discussion of 
the implication of this line of work. 

PACS numbers: 47.27.Gs, 47.27.Jv, 47.80.Fg 


I. INTRODUCTION 


Flowing fluids are ubiquitous in many natural and indus¬ 
trial situations. The running of water out of the faucet in 
our kitchen or the intimidating roar of a destructive hurricane 
provide examples involving the two most abundant fluids on 
earth: water and air. From elementary physical principles, the 
description of fluid motions is based on the competition be¬ 
tween the inertia of fluid particles, and the diffusion of mo¬ 
mentum by viscosity. As a result, the physical properties 
of the flow are characterized by the dimensionless Reynolds 
number Re = ULjv, where L and U are the typical length 
and velocity scales of the flow and v is the kinematic viscosity 
of the fluid. The Reynolds number can be regarded as the ratio 
of the viscous time scale L'^/v and the flow time scale LjU 
and therefore measures the relative importance of the inertial 
effect, which tends to drive the flow to become unstable, and 
the damping by the viscosity. For a flow at small Reynolds 
number, viscous diffusion is fast, so viscosity damps out flow 
disturbances. This situation is referred to as “laminar”. In 
such flows, energy dissipation transfers mechanical energy 
into heat, which can be readily seen from the few available ex¬ 
act solutions of the Navier-Stokes equations Q] . On the other 
hand, when the Reynolds number is large, the inertial effects 
dominate and the flow appears to be much more irregular, even 
in the absence of any externally imposed time-dependence. 
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Such flows are called “turbulent”. In three dimensional sit¬ 
uations, turbulent flows are “rough”, in the sense that they 
develop strong variations of the velocity field over very small 
scales, or equivalently, very large velocity gradients. In these 
regions, the viscosity is important. The irregular nature of tur¬ 
bulent flows leads to a fast and seemingly erratic motion of 
small particles transported by turbulence. The work reviewed 
here shows that the fundamental properties of the flow, such as 
the irreversibility induced by the energy dissipation, manifest 
themselves in the motion of small particles. 

Most macroscopic flows in nature and technology are tur¬ 
bulent. This is a consequence of the very small values of the 
viscosities of the most common fluids, such as air and wa¬ 
ter, which leads to large Reynolds numbers, even at modest 
length and velocity scales. For example, an adult walking at 
a moderate pace creates an air flow with a Reynolds number 
of approximately 5 x 10^ around him/her, and the flow of tap 
water in our kitchen can easily reach a Reynolds number of 
10^, both of which are well in the turbulent flow regime. 

Because of their rapid erratic motion, turbulent flows 
strongly enhance mixing. This we know well from stirring 
water in order to dissolve sugar added in it. The same princi¬ 
ple we apply when we rapidly mix fuel and air in combustion 
engines. This of course comes at the expense of energy: The 
intense, erratic turbulent flow needs to be maintained by exter¬ 
nal driving or pumping. Sometimes one would like to avoid 
turbulence. Examples can be found in fluid drag on trains, 
cars, ships and airplanes, or in pipe flows, where a much larger 
pressure drop is required to pump a turbulent flow through a 
pipe than a laminar flow at the same flow rate. Thus a better 
understanding of turbulence and turbulent flows could allow 
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us, on the one hand, to mix fluids more efficiently and, on the 
other hand, to reduce the drag in technical applications. 

From the point of view of fundamental physics, turbulence 
is an emblematic example of non-equilibrium systems, whose 
description is notoriously challenging. In particular, novel 
concepts to master the underlying complexity of turbulent 
flows are yet to be discovered. These concepts could then 
also provide deeper insights into other non-equilibrium prob¬ 
lems. One of the unifying concepts to address these problems 
is the breaking of detailed balance. Irreversibility in turbu¬ 
lence implies that the transition probabilities from a state A 
to a state B, and from the state B to the state A are not equal, 
contrary to what happens in equilibrium systems in statistical 
mechanics M. 

Like other out-of-equilibrium systems, turbulent fluid flows 
are irreversible. The kinetic energy of the flow is always dis¬ 
sipated and a constant supply of kinetic energy is necessary 
to maintain a turbulent flow. Whereas it is a simple mat¬ 
ter to understand dissipation (hence irreversibility) in laminar 
flows, it is much more challenging to identify the features of 
the turbulent flow motions that reveal irreversibility. In turbu¬ 
lent flows, the nonlinearity plays a crucial role, and leads to 
chaotic motion, with subtle statistical properties. A specific 
property of turbulent flows is that the scales at which energy 
is supplied, either from external forcing or from flow instabil¬ 
ities, are vastly different from the scales at which the energy 
is dissipated. In three-dimensional (3D) flows, the energy in¬ 
jection is at large scales while the viscous dissipation domi¬ 
nates at small-scales iUHl. In two-dimensional (2D) flows, 
the energy is supplied at small scales and then dissipated by 
friction at large scales The scale-separation between 

the energy injection and energy dissipation implies that there 
is, on average, a flux of kinetic energy through spatial scales 
between the energy injection scale and the energy dissipation 
scale, i.e., in the so-called inertial range, and the direction of 
this average energy flux cannot be reversed. This is the under¬ 
lying reason why turbulent flows are irreversible. There are 
fundamental differences between fluctuations in steady-state 
turbulence and fluctuations in equilibrium systems ifl^ . 

Remarkably, the celebrated Kolmogorov’s 4/5-law shows 
that in statistically stationary turbulent flows, the energy flux 
can be measured from a single snapshot of the velocity 
field ifTlI [l^ , without any explicit reference to its temporal 
evolution. The 4/5-law can be generalized to both 2D and 3D 
cases as (see e.g. Ref. ifT^ ) 
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where u(x, t) is the velocity field, r is the separation vector 
whose magnitude is r, d is the spatial dimension, and e is 
the turbulent energy flux, defined as a signed quantity, with 
the convention that e > 0 when energy is transferred towards 
small scales and e < 0 for energy flux toward large scales. 
The symbol (•) refers to an ensemble average, obtained by 
averaging over all possible flow realizations. For statistically 
homogeneous flows, it is the same as averaging over time t 
(under the assumption of ergodicity), or over space x. Equa¬ 
tion ([T]) establishes that the third moment of the longitudinal 


velocity increments, SrU = (^u(x -t- r,t) — u(x,/)j • (r/r), 
differs from zero, which is a consequence of the existence of a 
flux of energy through scales. This equation is valid when the 
distance r between the two measured velocities is in the iner¬ 
tial range, i.e., when is r i^, where is and iz, are the 
smallest and the largest length scales of the flow whose phys¬ 
ical meanings depend on the spatial dimension of the flow. 
For 3D turbulence, ig = rj, the well-known Kolmogorov scale 
given by the balance between the viscous dissipation and the 
energy flux e (remember that e can be either positive or neg¬ 
ative, depending on the nature of the transfer of energy), and 
iL is the forcing scale. In the 3D turbulence case, e > 0 and 
the flux is from the large to the small scales and it is called 
direct energy cascade. On the other hand, for 2D turbulence, 
the energy flux is negative: e < 0 and it is called inverse en¬ 
ergy cascade because the kinetic energy is from the forcing 
scale if, which is the small scale ig, to the large scale that 
is determined by the balance between large-scale friction and 
the energy flux. 

We stress that in Eq. O, both u(x 4- r) and u(x) are mea¬ 
sured at the same time t, i.e., in principle using only a sin¬ 
gle “snapshot” of the velocity field. This property leads to 
the following paradoxical situation. Consider a sequence of 
velocity fields, which correspond to a solution of the Navier- 
Stokes equation. According to Eq. [l] the precise order of the 
snapshots is immaterial in the determination of the rate of en¬ 
ergy dissipation. In particular, simply reversing the order of 
the sequence, i.e., changing t ^ —t does not seem to affect 
the determination of the energy dissipation rate, which is ob¬ 
tained from the individual velocity field (Eulerian statistics) 
alone. Thus, we can measure the energy flux correctly but 
cannot detect the time-irreversibility, although the energy flux 
is the cause of the irreversibility. 

To understand this paradoxical situation, we investigate 
here the motion of fluid particles moving with the flow, i.e., 
Lagrangian statistics. In the last two decades, significant 
advances in measurement techniques and computing tech¬ 
nologies made it possible to obtain well-resolved Lagrangian 
statistics at high Reynolds numbers in both experiments |fl4| - 
1^ and numerical simulations Unexpectedly, as ex¬ 

plained in Section UlI A1 the Lagrangian structure functions 
turn out to be completely insensitive to the fundamental irre¬ 
versibility of the flow 112^ . thus calling for new concepts and 
ideas. In the following, we review briefly these recent results, 
with the focus on the relation to the irreversibility of turbu¬ 
lence. 


II. MULTI-PARTICLE STATISTICS: THE ROLE OE 
ENERGY ELUX 


We start with multi-particle Lagrangian statistics, i.e., by 
following more than one fluid particle in a turbulent flow. To 
this end, we consider elementary sets of particles, and study 
the flow perceived by the particles, and its dependence on the 
characteristic distance between the particles. We show that 
when the inter-particle separations are in the inertial range, the 
energy flux through scales results in a measurable difference 
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between the statistical properties of the motion forward and 
backward in time. 


A. Relative dispersion and energy considerations 


In a turbulent flow, two fluid particles are, on average, mov¬ 
ing away from each other. This can be quantified by mea¬ 
suring the mean squared stance between the two particles, 
{R‘^{t)) = (|Xi(f) - X 2 (f)P)flo where Xi(f) and X 2 (f) are 
the positions of particles 1 and 2 at time t. The subscript Rq in 
the definition of the mean squared distance between two par¬ 
ticles, refers to the imposed condition that at f = 0 

the initial distance i?(0) = |Xi(0) — X2(0)| is equal to Rq. 

How changes with time quantifies the relative dis¬ 

persion of a particle pair. The problem has been extensively 
studied since the pioneering work by Richardson llz^ . who 
observed that in turbulent flows, in a meteorological con¬ 
text, the mean squared distance between two particles grows 
with time as i.e., (i?^(f)) oc Paradoxically, although 
Richardson law is easy to justify theoretically, it has proven 
very difficult to observe this regime of turbulent dispersion in 
any other well-controlled laboratory flows, or in direct numer¬ 
ical simulations, and r nany investigations have been devoted 
to this phenomena lfl4lll8ll2^l28l - l^ . 

In general, the separation (i?^(f)) can be expressed as 

(i?2(f)) =Rl + 2 [ [ {6u{t') ■ du{t"))dt"dt', (2) 

Jo Jo 

where Su{t) = Ui(f) — U 2 (f) is the relative velocity be¬ 
tween the two fluid particles. Equation (|2]l thus relates the 
relative dispersion of two particles to the Lagrangian correla¬ 
tion of relative velocity, (5u(0)(5u(r)). In fact, considering 
{(5u(0)5u(t)) in the limits of small values of r (r —>■ 0), and 
T large compared to the velocity correlation time (r oo), 
leads to interesting information on the time-dependence of 
(i?^(f)). When T —>■ cxD, the relative velocity (5u(r) becomes 
independent of its initial value (5u(0), so the double integral 
in Eq. (|2]l is linear in t. This corresponds to a diffusion-like 
regime that is similar to the turbulent diffusion of single ar¬ 
ticles first discussed in the seminal work of G. I. Taylor 11^ . 
In the opposite limit of very short times, r —0, we can ex¬ 
pand the integrand in Eq. (|2]i at t = 0 in power series of t and 
integrate to obtain; 

(R^(t)) =R^ + (Su^(0)}t^ + (Su(0)-Sa(0)}t^ + 0(t^), (3) 


where Su^(O) is a simplified notation for (5u(0) • (5u(0) and 
da(t) = ai (t) — a2 (t) is the relative acceleration between 
the two particles. Equation Q shows that as long as the 
initial separation Rq is non-zero, the initial velocity differ¬ 
ence does not vanish and is dominated by the 

term at very small times. The Richardson regime, there¬ 
fore, can only exist for some intermediate time t after mark- 

Eurthermore, we note that 

is the rate of change of the ki¬ 


ing the two particles 

Su{t) ■ 5a{t) = ^ \5v?{t) 


netic energy in the relative motion between the two particles. 


It has been shown that for separation Rq in the inertial range, 
this rate of kinetic energy change is related to the turbulent 
energy cascade through scales llo[l4l^ : 

= = (4) 

in which the averaging is taking over all particle pairs sepa¬ 
rated by a distance Rq in the inertial range. As we mentioned 
before, the energy flux e is positive for 3D flows. Hence 
in 3D turbulence the kinetic energy in the relative motion 
conditioned on a given separation Rq between particles ini¬ 
tially decreases. It increases at later times, consistent with 
the faster than increase of {R^{t)) at later times. This un¬ 
expected consequence of the energy cascade has been con¬ 
firmed, both in numerical simulations and in Lagrangian 
particle tracking experiments s. Note that substituting 
Eq. dUl into Eq. Q gives a negative term for {R^ (f)), which 
should not be confused with the positive coefficient in the 
Richardson dispersion law expected at later times. 

In the problem of mixing of a passive scalar, a proper mod¬ 
eling of the fluctuations of concentration rests on understand¬ 
ing how two fluid particles arrive at a given distance apart, 
or in other words, how {R?{t)) changes when t < 0. This 
amounts to tracking the motion of particles backward in time. 
For backward dispersion, it is also expected that a Richardson- 
like regime exists, i.e., f)) oc for intermediate time 

t, but with a larger coefficient, which means that backward 
dispersion is faster than forward dispersion ll4d1 - l43l] . An in¬ 
teresting observation is that for incompressible Navier-Stokes 
turbulence, both Eq. (|2|l and Eq. Q are also valid for f < 0. 
We note that in 3D turbulence, the term in Eq. Q is positive 
for f < 0, which implies that even for short times, backward 
dispersion is faster than forward dispersion. This is a conse¬ 
quence of the fact that the kinetic energy of the relative mo¬ 
tion of particle pairs, followed backwards in time, increases at 
short times, contrary to what happens when following particle 
motion forward in time. 

This property is a manifestation of time irreversibility of 
turbulence. In principle, it can be used to detect the “arrow- 
of-time”, while following many particle trajectories in a tur¬ 
bulent flow llril] . This manifestation of time-irreversibility in 
the relative dispersion between two fluid particles, d{R?)/dt, 
ultimately rests on the relation between the turbulent energy 
cascade and the rate of energy change in the relative motion 
expressed by Eq. dUl. This means that if we know the deriva¬ 
tive of the kinetic energy in the relative motion, d{^5v?)/dt, 
or the velocity held and its time derivative dSu/dt, then we 
can also detect the “arrow of time”. On the same basis, other 
Lagrangian quantities that combine both relative velocity and 
separation R can also be formed with the property that their 
time derivatives are sensitive to whether the “arrow of time” 
is flipped or not S. A further observation is that for other 
flows that do not satisfy incompressible Navier-Stokes equa¬ 
tions, Eq. (01 might not remain the same for t > 0 and f < 0, 
i.e., there could be an anomaly in Lagrangian velocity statis¬ 
tics. For example, it has been shown that for the compressible 
Burgers equation, taking d{Su^)/dt at f = 0 from the t < 0 
and t > 0 side give different values, which is due to the for- 
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mation of shocks in Burgers turbulence when time is running 
forward This is clearly a stronger manifestation of time 
irreversibility. 

While Eq. dUi has been derived for both 2D and 3D tur¬ 
bulence, its validity has so far been verified numerically and 
experimentally only for 3D turbulence, but not for 2D turbu¬ 
lence. The physics of energy cascade is completely different 
in 2D turbulence compared to that in 3D turbulence. 

In particular, the energy flux in 2D is towards larger scales, 
and consequently the kinetic energy in the relative motion be¬ 
tween fluid particle pairs is expected to increase initially. Con¬ 
fronting this prediction with numerical and experimental data 
could be an interesting and important work for the future. 

B. Shape deformation and structure of the flow 

While the previous section was devoted to the relative mo¬ 
tion between two particles, we now turn to the Lagrangian 
statistics involving more than two particles. The description 
of a set of points requires not only a size, such as the mean dis¬ 
tance between the particles, but also extra variables describing 
the shape of the set of points. The shape evolution provides in¬ 
teresting information on the local (topological) structure of the 
turbulent flow, which cannot be obtained from the study of the 
mean separation between pairs of particles alone. To explore 
the flow topology, one needs to follow at least 3 particles in a 
2D flow and 4 particles in a 3D flow. It has been observed that 
the evolution of initially isotropic objects (equilateral triangles 
or regular tetrahedra) in turbulent flows differs from that in a 
Gaussian velocity held. Qualitatively, the shapes obtained in 
a turbulent flow are more elongated at intermediate times than 
expected by using a flow with Gaussian statistics, or before all 
particles are widely sraarated so their velocities become inde¬ 
pendent 1^1^ 1^1^ . To see how this is related to flow 
topology, one can define an effective local velocity gradient 
M perceived by the set of particles 1^1^14^ : 

M = g-iW, (5) 

where the matrices g and W are defined as 

N 

9ij = 

and 

N 

a—1 

where N is the total number of particles in the set used to 
define the perceived velocity gradient M, 

1 ^ 

x'“ = x“-V x“ (8) 

N ^ ’ 

Oi—\ 

and 

1 ^ 

u'“ = U“ - - ^ U“ (9) 

' a=l 


are the position and velocity of particle a relative to the cen¬ 
ter of the particle set. It is easy to show that the perceived 
velocity gradient M given by Eq. (|5]l is the least square fit of 
a linear velocity held from the velocities at the particle po¬ 
sitions. When the separations between the particles are very 
small, in the range where viscous effects dominate, the per¬ 
ceived velocity gradient M given by Eq. (|5]) recovers the true 
velocity gradient. When the particle separations are in the 
inertial range, M probes the inertial range dynamics, which 
is the main motivation to study M 1^ . Other effective lo¬ 
cal velocity gradients similar to M have been proposed and 
the information on flow topology obtained are also compara¬ 
ble iOtI Hi- 

To probe the flow topology, it is helpful to decompose 
the velocity gradient M as a sum of a symmetrical part, 
S = (M -t- M^)/2, which represents the rate of strain 
(local stretching or compression) of the flow, and an anti- 
symmetrical part, r2 = (M — M^)/2, which represents the 
local rotation (by construction, S -f $1 = M). This decom¬ 
position is in fact unique. The strain S and the rotation 
interact with each other, which forms the rich dynamics of 
turbulent flows. In particular, on average the local rotation 
rate is constantly amplified because of the action of the strain, 
a phenomenon called “vortex stretching”, which is eventually 
compensated by the viscous dissipation. Early studies of the 
true velocity gradient in turbulent flows have revealed that the 
stretching of the vorticity is closely tied to the statistics of the 
eigenvalues of the strain S. Namely, among the three eigenval¬ 
ues, which are all real because S is symmetric by definition, 
the intermediate eigenvalue is predominately positive ifsill^ . 
This has been verified in numerical simulations and experi¬ 
ments ifs^ I 54 I 1 and has stimulated further studies on the dy¬ 
namics of velocity gradients ifssl - l?^ . A natural expectation 
from vortex stretching is that the vorticity vector would be 
preferentially aligned with the direction corresponding to the 
largest eigenvalue of the strain, which represents the strongest 
stretching. On the other hand, numerical and experimental 
data show that at any given instant, vorticity is preferentially 
aligned with the intermediate eigenvalue of the strain, which 
corresponds to rather mild stretching if^ 1^ . This inter¬ 

esting observation has been studied extensively in subsequent 
research (see M for a detailed discussion). 

When the size of the particle cluster used to obtain M from 
Eq. 0 is larger than the viscous range, M differs from the 
true velocity gradient and provides a way to probe the flow 
property in the inertial range of scales. It has been observed 
that M obtained in this way shares qualitatively many proper¬ 
ties of the true velocity gradient. In particular, the intermedi¬ 
ate eigenvalue of the rate of strain S is predominately positive 
and instantaneously the vorticity is aligned with this interme¬ 
diate eigenvalue These properties of M, which 

are closely related to the inertial range dynamics and hence 
to energy cascade, are expected to lead to irreversibility in 
Lagrangian multi-particle statistics. Eor example, the short- 
time deformation of an initially isotropic tetrahedron formed 
by four fluid particles is governed by the eigenvalues of the 
perceived rate of strain. Therefore, the non-zero average of the 
intermediate eigenvalue of S implies that the shape evolution 
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of a tetrahedron differs when followed forward or backward in 
time . The Lagrangian view also shows that the perceived 
vorticity vector indeed tends to align with the largest eigen¬ 
value of S, in the sense that the vorticity vector turns to the 
initially strongest stretching direction, but with a time delay 
such that at any given instant the vorticity is observed to pref¬ 
erentially align with the intermediate eigenvalue ll^. This 
property is also found for the true velocity gradient "lilli. 
The observed alignment process of vorticity with the rate of 
strain at a given time will completely differ when following 
flow trajectories backward in time. We note that when the 
length i?o characteristic of the set of points is in the inertial 
range of scales, the dynamical processes are essentially self¬ 
similar. This can be seen by rescaling time by the time scale 
to — (i?g/e)^/^, which can be viewed as the eddy-turnover 
time at scale Rq ll^l4l.l^ . 

A different but related interesting question is how rigid par¬ 
ticles with given shapes see the turbulent flow. For small rod¬ 
like particles, it is surprising to observe that they align with 
local vorticity and hence preferentially with the intermedi¬ 
ate eigenvector of the rate of strain 1163146^ . New results are 
available concerning the coupling between translation and ro¬ 
tation of neutrally buoyant particles with other shapes, such as 
large spheres ellipsoids or other anisotropic 

shapes iill]. It would be very interesting to And how their 
dynamics are related to the irreversibility of the flow. 

In summary, we note that for multi-particle Lagrangian 
statistics, the distance between particles defines a natural 
length scale of the problem and the energy cascade process 
in turbulence inherently causes the observed irreversibility. 
Stochastic models have been widely used to describe the 
multi-particle dispersion process and many aspects of the ob¬ 
served multi-particle statistics can be recovered by these mod¬ 
els iFt^. It^ . As almost all these models have time-reversible 
dynamics built in, one should be cautious not to push these 
models beyond the range in which they are valid. 


III. SINGLE-PARTICLE STATISTICS 


In situations involving several particles, we can naturally 
introduce one length scale (or more) in the problem, therefore 
permitting to establish a relation with the Eulerian correlation 
functions of the velocity held, and hence the energy flux or 
energy dissipation. In contrast, following only one fluid par¬ 
ticle in a turbulent flow does not give rise to an unambiguous 
identification of a length scale. For this reason, new ideas 
and concepts are needed in order to understand the flow prop¬ 
erties from the statistics of single particle trajectories only. 
Recent progress provides new insights into these interesting 
questions. 


A. Velocity structure functions 

The most studied single-particle statistic quantities are the 
Lagrangian velocity structure functions, i.e., the moments of 


the velocity increments following a fluid particle: 


Snir) = ((<5rM)"') = {[uit + t) - u(f)]”), (10) 

where u{t) is one component of the particle velocity at time 
t along a direction, being understood that for homogeneous 
and isotropic turbulence, the choice of the component does 
not matter. By analogy with the Eulerian velocity increments, 
one may surmise that the Lagrangian velocity increments SrU 
depend on the turbulent energy dissipation rate e. Further¬ 
more, if the time lag r is much larger than the viscous time 
scale but smaller than the largest time scale of the flow, it is 
tempting to postulate that the statistics of the velocity incre¬ 
ments are universal and independent of viscosity. Simple di¬ 
mensional analysis then leads to the scaling 6rU ~ 
and hence «S'„(r) (er)”/^ (ll [32. [H, Hi ■ Available ex¬ 

perimental and numerical data show that the dependence of 
S'„(t) on r has very little to do with the expected scaling be¬ 
havior dill in. Various theories have been proposed to 
explain the observed deviations li^ with the multi¬ 

fractal model being the most popular 175, 7^ (see also a 
recent summary in Ref. lislll l. 

Among the Lagrangian structure functions, the second or¬ 
der, obtained by taking n = 2 in Eq. (fTOb , is of special interest 
because according to the dimensional argument it is propor¬ 
tional to the energy dissipation rate e itself, so the average in 
Eq. ([Tol l is not affected by the strong fluctuations in e, which is 
known to lead to corrections to scaling in the case of the spa¬ 
tial structure functions (intermittency corrections) 

Based on these considerations, the scaling S 2 (t) oc er is ex¬ 
pected to be exact, just as the linear scaling of the third-order 
Eulerian velocity structure function predicted by Eq. ([T]l (the 
4/5-law). This expectation is summarized by the following 
relation: 


S 2 (t) = ((Sruf) = Coer, (11) 

where Co is expected to be a universal constant of order 
unity 11^ 1^ \lm . The best available data, from state-of-the- 
art experiments and numerical simulations, however, does not 
support the scaling suggested by Eq. (fTTT i H^lllll. If any¬ 
thing, the values of Co observed for 3D turbulence are found 
to increase with the Reynolds number of the flow and are ap¬ 
proximately 7 for the highest Reynolds numbers measured so 
far lH^IHl]. In the case of 2D turbulence, where Eq. (fTTl l is 
also expected to hold, the values of Co increase much faster 
with the Reynolds number. The results of numerical simula¬ 
tions at the highest available resolution suggest values of the 
order ~ O(IO^), without any indication of saturation ||2^ . 
This casts a serious doubt on the scaling predicted by using 
dimensional arguments. 

In fact, the assumption that the statistics of velocity incre¬ 
ments following a fluid particle, 5rU, depend on e, the en¬ 
ergy flux through spatial scales, is questionable. That assump¬ 
tion is directly inspired from Kolmogorov’s hypotheses on the 
statistics of Eulerian velocity increments. As already noticed, 
establishing a connection between Lagrangian and Eulerian 
statistics requires the introduction of a length scale into the 
structure functions S'„(r), which is achieved by assuming that 
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t6tU is equivalent to the separation r in Eulerian statistics 
and that SrU scales as the Eulerian velocity difference S^u. 
Substituting 6rU ^ 6rU and r ^ rSrU into Eq. ([TJ leads to 
{{Sru)'^) ^ er, hence to a formal justification of Eq. (fTTl i. 

The formal analogy between the Eulerian velocity incre¬ 
ments, SrU, and the Lagrangian velocity increments SrU, 
through the use of r Ri rSrU should, however, be taken very 
carefully. The Eulerian statistics of SrU are mostly deter¬ 
mined by turbulent eddies of size r. When estimating the La¬ 
grangian increment St-u using the Eulerian increment SrU with 
r w rSrU, it should be noticed that the time t necessary for a 
particle to travel up to r, t ^ rjSrU ^ (Sru)'^/e, is in fact the 
life time of an eddy of size r. That is to say, for a fluid par¬ 
ticle to move a distance r, its velocities at the start, u{t), and 
at the end, u{t + r), are unlikely to be the result of the same 
eddy of size r. This essential dissimilarity between SrU and 
SrU leads to very different properties between Eulerian and 
Lagrangian statistics. More generally, this feature highlights 
statistically stationary turbulence, as an ultimate example of 
non-equilibrium steady state, far from equilibrium ifl^ . 

These, and other considerations led Ealkovich et al. 
to question the validity of Eq. They pointed out that 

the statistics of SrU, including all Lagrangian velocity struc¬ 
ture functions, are symmetric under the transformation of 
t ^ —t, therefore being unable to pick up the fundamental 
time-irreversibility of the flow. Therefore, there is no fun¬ 
damental reason to relate the statistics of SrU to the energy 
flux, which is the cause of the time-irreversibility of turbulent 
flows. 

In summary, there is strong motivation to consider other 
Lagrangian statistics that reveal the irreversible nature of tur¬ 
bulent flows. That is the topic we cover in the next subsection. 

B. Kinetic energy increments and instantaneous power 

An interesting recent discovery is that the change of kinetic 
energy following individual fluid particles can be used to de¬ 
tect the “arrow of time”. In Ref. 1841] . it was observed in ex¬ 
periments and numerical simulations that the third moments 
of the kinetic energy change, SrW = [u^(t) — u^(0)]/2: 

+ r)/2 - u\t)/2f) (12) 

are negative for time lags r positive, but smaller than the ve¬ 
locity correlation time (the largest time scale of the flow). This 
implies that the probability distribution of the instantaneous 
power, p = limT-_>.o SrWjr = u-a, is negatively skewed. The 
origin of this skewness can be traced back to the observed ten¬ 
dency of fluid particles to gain kinetic energy slowly, but lose 
it more suddenly. This provides a way to identify the arrow 
of time, as flipping t ^ —t would lead to the exact oppo¬ 
site; particles would gain energy faster than they dissipate it. 
The negative skewness of the distribution of p was observed 
for both 2D and 3D turbulence (at least for 2D turbulent flows 
that were agitated with forces short-correlated in time), i.e., 
independent of whether energy flows towards larger or smaller 
scales. Erom the more general point of view of energy ex¬ 
change, kinetic energy is dissipated into heat in an irreversible 


way in turbulent flows, both in 2D and 3D. In that sense, the 
qualitative similarity between the statistics of p and SrW in 
both 2D and 3D turbulences, once expressed in terms of the 
energy dissipation, e, may not be so surprising. Available data 
support that the moments of p/e follow, to a good approxima¬ 
tion, a power law dependence on the Reynolds number of the 
flow, with an exponent independent of the spatial dimension. 
This suggests that the third moment — (p^) /e^ can be used as 
a measure of irreversibility. Moreover, the scaling of the third 
moment — (p^)/ can be qualitatively explained by assuming 
that — (p^) /e^ is dominated by the extreme events of negative 
p with large magnitudes, i.e., events when fluid particles lose 
kinetic energy very rapidly, an argument pictorially alluding 
to “flight-and-crash” events ifsdl] . This skewed distribution of 
p also manifests itself in the negative skewness of the kinetic 
energy change associated with single velocity component in 
a 3D turbulent flow ll^ . and in the negative skewness of the 
longitudinal Lagrangian velocity increments 1^ . 

While the skewness of the instantaneous power p is nega¬ 
tive for both 2D and 3D turbulence, and in this sense, seems 
to be insensitive to the very different physical mechanisms of 
cascade in these two cases, one may nevertheless ask which 
quantity reflects the difference in the dynamics in 2D and 3D 
flows. To answer that question, one can decompose the in¬ 
stantaneous power p into 

p = u • a = —u • VP -f u • f + u • D, (13) 

where —VP, f, and D are the pressure gradient, external 
forces, and dissipative forces, respectively. In 3D flows, the 
dissipative forces consist of the viscous forces alone, D = 
I'V^u; while for 2D flows, the dissipative forces include both 
viscous forces and friction forces, D = — au, where 

a > 0 is the linear friction coefficient. Numerical simulation 
data show ifs^ that in both 2D and 3D flows, the magnitude 
of the pressure gradient term —u ■ VP overwhelms all other 
terms and determines the magnitude of p, but the contributions 
to the third moment of p are more subtle and show interesting 
differences between 2D and 3D flows. 

In 2D flows, the pressure gradient term is also negatively 
skewed, and it contributes to nearly 2/3 of (p^), with the 
other dominant contribution being provided by the correla¬ 
tion between the pressure gradient and the friction, {(—u ■ 
VP)^(—au • u)). In 3D flows, the situation is completely 
different; the skewness of the pressure gradient term is very 
small, even slightly positive, so its direct contribution to (p^) 
is very small and of opposite sign, compared to (p^). The 
dominant term that contributes to (p^) is the cross term be¬ 
tween the pressure gradient and the viscous forces, {(—u ■ 
VP)2(iyu-V^u)) 0. 

Therefore, the pressure gradient term acts very differently 
in 2D and 3D flows. In 2D, it behaves according to naive ex¬ 
pectation, insofar as it provides the dominant term of the fluc¬ 
tuations, and contributes significantly to the observed asym¬ 
metry of the distribution of p. In comparison, although the 
pressure term in 3D also provides the main contribution to the 
variance of p, it hardly provides any significant contribution 
to the observed asymmetry of the distribution of p. 

Note that in homogenous flows, the pressure gradient term 
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alone does not change the total energy in the flow since 
(u • VP) = 0, i.e., the pressure gradient term merely re¬ 
distributes kinetic energy within the flow. The different role 
played by the pressure gradient term would imply that the way 
of energy redistribution is different in 2D and 3D flows. In¬ 
deed, the averaged value of — u • VP conditioned on the ki¬ 
netic energy of the particles reveals that in 2D flows, parti¬ 
cles get as much energy from pressure gradient forces as they 
lose it, independently of their velocity: (—u • VP|u^) = 0 
for all u^. In contrast, in 3D flows, the mean pressure con¬ 
tribution conditioned on the velocity is negative for particles 
with small velocities: (—u • VP|u^) < 0 for < 2(u^), 
and it is positive for particles with large velocities. This im¬ 
plies that the pressure gradient term takes kinetic energy away 
from slow particles and gives it to fast particles isl . With¬ 
out other terms to stop this action, the pressure gradient term 
alone could potentially drive the flow into singularities. This 
observation might provide new insight into the long-standing 
Millennium problem on the regularity of the Navier-Stokes 
equations 1 ^1^ . 

The decomposition of the instantaneous power p in the form 
of Eq. (foi l is certainly not unique. A possible alternative con¬ 
sists in decomposing the fluid acceleration into a local part 
sll = du/dt and a convective part slc = u-Vu. This decom¬ 
position separates the effect of the flow seen by the particle as 
resulting from the time variation of the velocity field locally, 
a^, and from the advection by a time-independent (frozen) 
flow, ac. It has been noticed that the two components a/, and 
Slc cancel each other to a large extent 1^ . How their 
contributions to the instantaneous power p behave and what 
they reveal about the irreversibility of the flow are interesting 
problems for future study. 

The results discussed above are all for incompressible tur¬ 
bulence. As predicted by the study of the two-particle statis¬ 
tics governed by Burgers equation 1461] . in compressible flows, 
the irreversibility is expected to manifest itself in a stronger 
way. It would be interesting to confirm these predictions, by 
either experimental or numerical simulation data. 

Another direction worth exploring concerns the effect of 
particle inertia on the effects discussed here for fluid parti¬ 
cles. Particles whose densities differ from the fluid density or 
whose sizes are larger than the Kolmogorov scale do not fol¬ 
low the flow faithfully due to their inertia. There has been a 
wealth of literature on the dynamics of various inertia parti¬ 
cles ifl^ 193 - 1^ . It will be therefore interesting to carry out 
analyses similar to that in Refs, H^] to see how the irre¬ 
versibility of the flow is reflected in the dynamics of inertial 
particles. 

IV. SUMMARY AND DISCUSSION 

In this brief review, we discussed how Lagrangian statistics, 
obtained by following the motion of fluid particles in the flow, 
are sensitive to the intrinsic time irreversibility of turbulence. 
This irreversibility is a consequence of the property that, in 
turbulent flows, the kinetic energy is supplied into the fluid 
motion at a scale very different from the scale at which the en¬ 


ergy is dissipated. In high-Reynolds number flows, where the 
forcing scales and the dissipation scales are widely separated, 
energy is transferred at a constant rate through scales, over 
a wide range of scales. This energy cascade and dissipation 
process are irreversible and the issue is how these phenomena 
affect the motion of fluid particles. 

We emphasize that depending on the spatial dimension, the 
energy is transferred either to small scales or to large scales, 
to be dissipated either by viscosity (in 3D), or by friction (in 
2D). This results in very different physical mechanisms. The 
investigation of multi-particle Lagrangian statistics, obtained 
by following sets of particles, allows the identification of at 
least one length scale. As a consequence, multi-particle La¬ 
grangian statistics naturally sense the energy flux and reflect 
it in the change of kinetic energy associated with the relative 
motion, the dynamics of the perceived velocity gradients and 
the shape deformation of isotropic objects. Single-particle 
statistics, on the other hand, does not permit such an unam¬ 
biguous identification of a length scale. Previous attempts to 
construct a length scale from single particle trajectories and 
hence to connect single-particle statistics with Eulerian statis¬ 
tics that depend on the energy flux are fundamentally ques¬ 
tionable, and lead to incorrect predictions. While the statistics 
of the usual Lagrangian structure function do not permit to 
distinguish the arrow of time, the statistics of other quantities, 
such as the energy change following a fluid particle, do re¬ 
veal the irreversibility of the flow. In this case, irreversibility 
is due to the energy dissipation, i.e., kinetic energy is eventu¬ 
ally converted to thermal energy. Those statistics, therefore, 
show similar behavior in both 2D and 3D, despite the oppo¬ 
site directions of the energy flux in flows with different spatial 
dimensions. 

Most of the work reviewed here originated from obser¬ 
vations in physical experiments and numerical simulations 
that built on rapid progresses in experimental techniques and 
computational power. As more and more high-quality, high- 
resolution data are generated by the research community and 
are being shared with the whole community (e.g., in large pub¬ 
lic databases we expect that other properties of La¬ 

grangian statistics in turbulent flows will be discovered and 
their connection with the time irreversibility of the flow, or 
more generally, with the dynamical properties of the flow, be 
understood. Lrom a broader theoretical point of view, fluid 
turbulence is a well-known example of out-of-equilibrium 
system. How to relate what we learned from studying fluid 
turbulence to other non-equilibrium systems is another open 
area for future investigation. 
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